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BOUNDARY LAYER SOLUTIONS OF CHARGE CONSERVING 
POISSON-BOLTZMANN EQUATIONS: ONE-DIMENSIONAL CASE* * * § 

CHIUN-CHANG LEEt, HIJIN LEE*, YUNKYONG HYON®, TAI-GHIA LIN’, AND CHUN LIUH 


Abstract. For multispecies ions, we study boundary layer solutions of charge conserving Poisson- 
Boltzmann (CCPB) equations [50] (with a small parameter e) over a finite one-dimensional (ID) spatial 
domain, subjected to Robin type boundary conditions with variable coefficients. Hereafter, ID bound¬ 
ary layer solutions mean that as e approaches zero, the profiles of solutions form boundary layers near 
boundary points and become flat in the interior domain. These solutions are related to electric double 
layers with many applications in biology and physics. We rigorously prove the asymptotic behaviors of 
ID boundary layer solutions at interior and boundary points. The asymptotic limits of the solution val¬ 
ues (electric potentials) at interior and boundary points with a potential gap (related to zeta potential) 
are uniquely determined by explicit nonlinear formulas (cannot be found in classical Poisson-Boltzmann 
equations) which are solvable by numerical computations. 

Key words, charge conserving Poisson-Boltzmann equations, boundary layer, multispecies ions 
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1. Introduction Almost all biological activities involve transport in ionic solu¬ 
tions, which involves various couplings and interactions of multiple species of ions. Many 
complicated types of electrolytes involved in biological processes, such as those in ion 
channel proteins, certain amino acids (movable side chain) are crucial to the functions 
of these ion channels. The electrostatic properties involving multispecies (at least three 
species) ions can be fundamentally different to those with only one or two species [U |33] . 
To see such difference, we study charge conserving Poisson-Boltzmann (CCPB) equation 
for multispecies ions which is derived from steady state Poisson-Nernst-Planck systems 
with charge conservation law, and is the surface potential model for the generation of 
a surface charge density layer related to electric double layers (SO] [50] . For simplicity 
of analysis, we consider a physical domain xG (—1,1) with the simplest geometry, and 
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represent CCPB equation as follows: 


( 1 . 1 ) 


-e^6" = 


N 

E- 


Zieonii 




for a;G(—1,1), 


where is the total concentration of species i with the valence Zi, (j) is the (electri¬ 
cal) potential, Cq is the elementary charge, ks is the Boltzmann constant, and T is 
the absolute temperature. The parameter e= {e^Ur/{<PeS)) >0, where cq is the 
dielectric constant of the electrolyte, Ut is the thermal voltage, d is the length of the 
domain (—1,1), and S is the appropriate concentration scale (cf. [H]). Furthermore, 
ed is known as the Debye length and e is of order 10“^ for the physiological cases of 
interest (cf. [3)- Thus we may assume e as a small parameter tending to zero. Similar 
equations to (ED can also be obtained by the other variational method [53) . 

Under suitable scales on and e, we let —Oj’s be the valences of anions, i.e., 
ak = —Zk, k=l,- ■ ■ ,Ni and bj’s be the valences of cations, i.e., bi = zi, 1 = 1, - • •, V 2 . Then 
the total concentrations of anions and cations are approximately given as Ok ~ ruk 
(fc = l,...,7Vi) and Pi^mi (/= 1,-• •,V 2 ), respectively. Hence equation (11.11) can be 
transformed into 


Nl N2 

= y - y 

( 1 . 2 ) 

for xS(— 1 , 1 ), 

where Ofc’s and bi’s satisfy 1 < oi < 02 < • • • < ojvi and 1 < 61 < 62 < ■ • ■ < ■ 

Most of the physical and biological systems possess the charge neutrality (zero net 
charge). One may assume the pointwise charge neutrality i.e. at all points the anion 
and cation charges exactly cancel in order to make calculations easier in a free diffusion 
system (cf. [19] p. 319). Here we replace the pointwise charge neutrality by a weaker 
hypothesis called the global electroneutrality being represented as 

Afi N2 

(1.3) Global Electro-neutrality: E akO!k = y^bil3i, 

k=i 1=1 


which means that the total charges of anions and cations are equal, where —akS and 
bi's are the valences, and a^’s and PiS are the concentrations of anions and cations, 
respectively. Consequently, the CCPB equation EU may satisfy ESJ. 

When one deals with more general (realistic) situations, such as when there are 
more than two species involved in the solution, situations become more subtle and 
complicated. Note that the equation m has nonlocal dependence on (j)^ which is 
essentially different from the classical Poisson-Boltzmann (PB) equation as follows: 


^1 ^2 A « 

(1.4) e2^"(a;) = y for a:e(-l,l). 

Here ^’s and ^’s are bulk concentration of anions and cations, respectively. In equa¬ 
tion El, Q!fc’s and /3/’s are for total concentration of anions and cations, respectively. 
For notation convenience, we use the same notations afc’s and /3;’s in equations El 
and El, but with different physical meaning. In this paper, we shall show different 
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asymptotic behaviors of the CCPB equation (11.211 and the PB equation p.4ll for various 
constants Ni,N2,ak,oik,bi,/3i satisfying dn. The main goal of this paper is to compare 
the CCPB equation (11.21) and the PB equation (11.41) under the hypothesis (11.31) . Such 
a difference can be clarified in Theorems o and 11.31 see also, Remark ITtI 

Boundary effects are important in a wide range of applications and provide 
formidable challenges [521 Hi]. For CCPB equations, the main issue is how bound¬ 
ary conditions effect the solution values (electric potentials) at interior and boundary 
points. One may use the Neumann boundary condition for a given surface charge dis¬ 
tribution and the Dirichlet boundary condition for a given surface potential (cf. [T]). 
Here we consider a Robin boundary condition [iisaisoiiiaiiiiiTiiiiisT] for the 
electrostatic potential (() at x = ±l is given by 

(1.5) <(),(l)-b77e<()'(l) = <()J, <()<: (-!)-??£<()'(-!) = </>(( • 

where are extrachannel electrostatic potentials and r]^>0 is the coefficient de¬ 

pending on the dielectric constant [21137], and related to the surface capacitance. The 
parameter ratio ? 7 e = es/C's can be viewed as a measure of the Stern layer thickness, 
where es and Cs are the effective permittivity and the capacitance of the Stern layer, 
respectively (cf. [ 6 ]). Thus we may regard ^ as the ratio of the Stern-layer width to 
the Debye screening length. Similar discussion can also be found in m and m- To see 
the influence of ^ on the asymptotic behavior of (p^s^ we consider the limit limej^g ^ 
to be either a non-negative constant 7 or infinity. 



Fig. 1.1. Schematic picture of Robin boundary condition, <pc^Ve{4>i)x = (p^ atx = ±l, and 
the limit values t^lime^opeil), c = lime_>o<?(>e(®), ®€(—1,1) and (p = t — c. 


Suppose limej^o ^ = 7 he- ^ where 7 is a non-negative constant. Then we show 
that the solution of dOD with (11.51) satisfies lime 4 ,o</>e(±l) = and lime 4 ,o^£(a^) = c 
for x€ (—1,1), where c and t can be uniquely determined by (I1.16|) - (I1.18I) which imply 
that the value c is changed with respect to t (see Figure [TTT]) . Moreover, the potential 
difference C,=t — c is decreasing to 7 (cf. Theorem 11.31) . Note that as the parameter e 
goes to zero, the solution (p^ has a boundary layer producing the potential gap C,=t — c 
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affected by Stern and Debye (diffuse) layers and related to zeta potential (cf. [22] ) which 
plays an important role in ionic fluids. However, for the PB equation (II.41) . the value c 
must be zero which is independent of t and 7 (cf. Theorem ll.il) . This shows the difference 
of the CCPB equation dOj) and the PB equation (El which can also be observed 
by numerical experiments (See Figure |5T] and Table ISTl in Section |5|). Furthermore, 
numerical computations give several conditions to let the profile of function c to 7 
become monotone decreasing and increasing fFigure [5.21 and 15.31 in Section |5|) and non¬ 
monotone (Figure inn in Section |S]). 

In [3D], we studied the CCPB equation (11.21) for case of Vi = V 2 = 1, (ai,6i) = (1,1) 
and (ai,/3i) = (a,/3), i.e., the case of one anion and one cation species with monovalence. 
In this case, equation El can be rewritten as 


( 1 . 6 ) 

(1.7) 


e^(p'^{x) = ne{x)—pe{x) for a:S(— 1 , 1 ), 


n^{x) = 


o4'e{x) 




and p^(x) = 


/3e 






where ne{x) and Pe{x) represent (pointwise) concentrations of anion and cation species, 
respectively. When a = (3 holds (the electroneutral case), we had shown previously that 
lime 4 ,o?^e(a^) = lime 4 ,o 7 'e(a:) = f for a;G(—1,1). Moreover, the CCPB equation (ll.6l) - (ll.7l) 
and the conventional PB equation t^w”{x) = ^ have same asymptotic 

behavior (cf. Theorem 1.4 of [3D]). In order for the readers to compare those with 
the results in the current paper, most results of m are summarized in Appendix. To 
certain degrees, it also justifies why in many situations, PB equation provides more or 
less expected solutions. On the other hand, we consider the non-electroneutral case, i.e. 
a^/3. Without loss of generality, we assume a<f3 i.e. j_^ne{x)dx < f_^pe(x)dx which 
means that the total concentration of anion species is less than that of cation species. 
Then we prove that limej,o’^e( 2 ;) = lime^o 7 'e( 2 ;) = § for a;G(—1,1), but lime 4 ,one(±l) = 

0 < limgj,oe^Pe(=bl) = (cf- (11-251) '). This shows that electroneutrality holds true 

in the interior of (— 1 , 1 ), but non-electroneutrality occurs at the boundary points ± 1 . 
Furthermore, the extra charges are accumulated near the boundary points ±1 (see 
Theorem [T3D . 

The mixture of monovalent and divalent ions such as Na+, K+, Cl“ and Ca^+ plays 

the most important roles for vital biological processes. For instance, opening and closing 

of ionic channels is accomplished by escape or entry of Ca^+ into the channels (cf. [T 8 ]l. 

The voltage may depend on [Ca^+] the concentration of Ca^’*' (cf. [H]). Differences in 

ionic concentrations create a potential gap across the cell membrane that drives ionic 

currents (cf. [3D] P. 34). To see how the voltage i.e. (electrical) potential depends on 

[Ca^+], we may use the equation (11.21) with Vi = l, ^2 = 2, ai=bi = l and 62 = 2 to 

describe the mixture of Na"*" (or K'*'), Cl“ and Ca^"*" ions, where ai ^ [Cl“], /3i [Na"*"] 

and P 2 [Ca^+j. In Theorem 11.31 (ii), we prove that when the electro-neutrality holds, 

that is, ai =/?i-|- 2 /? 2 , the solution of EU satisfies \im(j)^{x) = c for a;G(— 1 , 1 ) and 

€—>■0 

cS (c*, 0 ) is uniquely determined by (11.161) and 

l-e3'^cosht _/3i _ [Na+] 
e'^sinhc /32 [Ca^+] 

where t = lim(^g(l) > 0, and c* = Mogsecht is a negative constant (see Remark ll.2l) . The 
6—^0 

formula (11.81) shows that the interior potential (voltage) c is increased if the boundary 
potential t is fixed and the ratio [Na+]/[Ca^+] is increased e.g. [Ca^+] is decreased and 
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[Na+] is fixed. Furthermore, Theorem 11.31 is also applicable to the other cases with 
multi-species ions including multivalent and polyvalent ions so the formula (II. 8|) can be 
generalized to 


(1.9) 


„ „(l-l-z'lc/- sinh(zt) ^ 
^ ^1 sinht ' 


2e‘^sinhc 


/32’ 


for oi = &i = 1 and b 2 = z>2 (see Remark 11.21) . Note that (|1.9I) shows how the value c 
depends on the value t. Such a result cannot be found in the PB equation (11.41) . 

1.1. Asymptotic behavior of the PB equation (|1.4l) - (|1.5ft The PB equation 
m with the boundary condition m can be regarded as the Euler-Lagrange equation 
of the energy functional 

(1.10) Ef^[u] = l + /(z,))dx+^[(4--«(-l))2 + (4+-«(l))^], 

for mG iJ^((—1,1)), where 

Ni N2 

(1.11) /(s) = ^afce“'''*-|-^/3ie“^''* for sGM. 

fc=i 1=1 


For the PB equation (11.41) with the boundary condition (jl.5|) . we study the asymp¬ 
totic behavior of the solution of m as e approaches zero. The boundary con¬ 
dition (ILSI) plays a crucial role on the monotonicity of 4><i- Here we consider three 
cases for the signs of (/)q and (f)Q-. (a) min{(()g ,(()g } > 0, (b) max {(/)q , (()q } < 0 and 

(c) min{^(j',(/)(C}<0<max{^Q ,0g }. Then the corresponding results are stated as fol¬ 
lows: 


Theoreivi 1 . 1 . Assume Y^^LiakO-k = Y^f2i^iPi- .bet G1,1))(~)C'^([—1,1]) be 

the solution of equation |1.^[ ) with the boundary condition 11.51) . Then 

(i) For ccG (—1,1), |<(>e(a::)| exponentially converges to zero as e goes to zero; 

(ii) //min{(()(j',(/)(]'}> 0, then (fe is convex on [—1,1] and 0 < (^^(a:) < max{(/)(]",(()|C} 
for a; G [— 1,1] . Moreover, there exists e* > 0 such that for 0 < e < e*, 4>e attains 
the minimum at an interior point o/ (— 1,1). 

(Hi) If ma,x{(j)Q ,(/)((} <0, then (fe is concave on [—1,1] and min{ </>][",(/)]]■}< (/)£(a:) < 0 
for x€[—l,l]. Moreover, there exists e*>0 such that for 0<e<e*, </>£ attains 
the maximum at an interior point of (—1,1). 

(iv) If min{(/)j]",(()|C} < 0 < max{(()Q',(()(]'}, then (/>£ is monotone on [—1,1] and 

min{ (()o , ((>0 } < < max{ 

(v) Ifliui — =7 and 0<7<cx), then = t uniquely determined by 

e 4 .o e e 4,0 

(1.12) \(fo -i]=' 7 {fit)-f{0)y/^ and min{0,(()(|'} <T< max{0, </)([■}, 

where f is defined by HI. 11}) . Moreover, t = t{j) is decreasing in"; if ((q >0 and 
increasing in "j if (j)(f <Q. 

1.2. The main results In this section we present the main results, which are 
about the asymptotic behavior of the solution of (O and (II3| as e goes to zero, in 
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our research of CCPB equation. The CCPB equation (11.211 with the boundary condi¬ 
tion (ESI) can be regarded as the Euler-Lagrange equation of the energy functional 


EM = 


(1.13) 


Cf' 


Ni 

i'\^dx + E afelog 

k=l 

^2 


N2 


^da; + ^/3dog 


-bill 


dx 


\21 




for 1,1)). The existence and uniqueness for the solution of (11.21) and (11.51) is 

the following proposition: 


Proposition 1.2. There exists a unique solution 1,1)) D 1,1]) of the 

equation ED with the boundary condition ED- 


The proof of the above Proposition 11.21 can be easily obtained from the arguments of 
[50] (see Appendix therein) and [5T] . 

Suppose 4>q =(j)Q = A and ^2^=1 o,kCtk = Then Proposition ll.2l implies the 

solution of ED and ED must be trivial and (/>e = A. To study the nontrivial solution 
of (II.2p and (|1.5I) . it is sufficient to assume pQ Replacing by pe + C for any 

constant C, one may remark that the equation ED is invariant. Consequently, without 
loss of generality, we may assume —pp = Pq > 0 hereafter. 

When Y^kLiO-kOik = Y^i^ibiPi, i.e., the global electroneutral case. Theorem 12.11 
shows that max 3 ;g[_i_i] |0(:(x)| is uniformly bounded to e and that p'^ exponentially 
approaches zero in (—1,1) as e tends to zero. Thus, it is expected that there exists a 
constant c such that all interior values of p^ tends to c as e goes to zero. Along with 
Lebesgue’s dominated convergence theorem, we have 

(1.14) lim / e“'''^'da; = 2e“''^ lim / = 2"'"^ 

407-1 407-1 

and then the energy functional (11.131) with u = pf^ approaches to the energy functional 
Ef^[pf\ as follows (up to a constant independent of pPj: 

(1.15) = f j {e^\p'f\^ + f{p^-c))dx+-^[{pp-pe{-l)f + {p-^-pe(l)f] , 

where / is defined by p.lll) . Here we have used lim^j^o 1^ =0 (by 

(|1.141) 1 and the approximation log(l -I- s) s with s= ^ J^^e°'’‘^'^^~'^^dx—1 to get 

logy iy e“*’^‘^'“‘^^da;-|-log(2e“'‘'^) — 1 as 0<e<Cl. 

Similarly, we have 

logy e“^'^'^'“°^da:-|-log(2e“^'‘^) — 1 as 0<e<Cl. 

Therefore, we show that in the case of global electroneutrality (II.3L the energy func¬ 
tional (11.131) approaches (11.151) . which has the same form as the PB energy func¬ 
tional (11.101) . 
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The asymptotic behavior of (j)tS at boundary a; = ±1 may depend on the scale of rj^. 

Here we study two cases for the scale of rj^^O: (i) limej^o^ = oo and (ii) limgj^o ^ = 7, 

where 7 is a nonnegative constant. Then the relation between the boundary value limits 

lim(/)e(±l) and the interior value limit c are demonstrated as follows: 

£^0 


Ni N2 

Theorem 1.3. Assume —(j)Q =(j>Q >0 and''^^akCtk = ''^^bil3i. Let (j>eGC°°{{—l,l))7\ 

_ fc=i 1=1 _ 

C'^([—1,1]) be the solution of equation il.tA} with the boundary condition i ff.5)) . Then 


lim())e(—1) = —t, lim0£(l)=t and \\m(j)f^{x) = c for a;G(—1,1), 
e4,0 e4.0 eJ'O 


where t and c are determined as follows: 

(i) If lim —= 00 , then c = t = Q. 

6 

(a) If lim —=7 andO<'y<oo, then (t,c) uniquely solves the following equations: 
elo e 

(1.16) </')l■-^ = 7 (/(^-c)-/( 0 ))^/^ 

(1.17) f{t-c) = f{-t-c), 

(1.18) |c|<t<.^+. 

Moreover, writing t = t{"f) and c = c{j) in (ii), we have 

(A) lim ^( 7 ) = lim 0 ( 7 ) = c* and lim t{'y)= lim 0 ( 7 ) = 0, where |c*| <</>]]" 

7—>-0 7—>-0 7^130 7—)-oo 

is uniquely determined by f{(j)Q—c*) = f{—(j)Q—c*). 

(B) <( 7 ) and ^( 7 ) — 0 ( 7 ) both are decreasing on (0,oo). 


Formally, using —>■ c in (—1,1) as e tends to zero, equation (II.2|) may approach to 

the following PB equation: 


(1.19) 


Ni 

e2^''(a;) = ^ 


1=1 ^ 


for a;G(—1,1), 


which may give results of Theorem 11.31 bv formal asymptotic analysis. However, in this 
paper, we focus on rigorous mathematical analysis and provide the proof of Theorem ll.3l 
in Section [2j 

Theorem o:i) shows that there is no boundary layer and —>■ 0 uniformly in 
[—1,1] as e 4,0 if lime Theorem ll.3H i]) assures the existence of boundary layers. 

Furthermore, Theorem o (ii-A) and (ii-B) represent the ratio of Stern screening 
length to the Debye screening length affects the boundary and interior potentials: (a) 
The decrease of 7 results in the increase of t — c (the potential difference between the 
boundary and interior); (b) If 7 —>- 00 , the potential difference t — c may approach zero. 
Notice that the formula ()1.12|1 is quite different from (II. 1611 - (II. 181) . This may show the 
difference between solutions of the CCPB equation (ll.2|l and the PB equation (ll.4|l . 


Remark 1.1. 

(a) Theorem 17.71 (ii) and (Hi) show that as (j)Q(j)Q>0, the solution of the PB 
equation ra may lose the monotonicity. However, the solution of the CCPB 
equation (m always keeps the monotonicity (see E.em,ark \2.2\ (i)). This pro¬ 
vides the difference between solutions of the CCPB equation and the PB 

equation O- 
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(b) For equation liL2\} . the values c (interior potential) and t (boundary potential) 
depend on each other and satisfy precise formulas However, 

for equation interior potential and boundary potentials (determined by 

il.l2]) ) are independent to each other. 

Remark 1.2. When Ni = \, N 2 = 2, ai = bi = l, 62 = 2 and ai = I3i+2j32, we may 

get U.8\) from U.16\) and Moreover, \1.9\) can also be derived from il.l6\) 

and for the case that Ni = l, N 2 = 2, ai = bi = l, b 2 = z>2 and ai=/ 3 i+z/ 32 - 

By dUgj and im, it is easy to check that ^ < 0 for t>0. Then c = c{t) can be regarded 
as an decreasing function to t>0. Consequently, by Theorem \1.3\ (iv), c is increasing 
to 7 . 


When 7Vi = 1, iV 2 = 2, ai = 61 = 1 and &2 = 2, further asymptotic behavior of (f>e near 
the boundary a; = ±l describing the boundary layers is stated as follows: 


Theorem 1.4. Assume Ni = l, N 2 = 2, ai=bi = l and 62 = 2 . Under the same hy¬ 
potheses of Theorem \2.1\ and Theorem \1.3Y ii). the asymptotic behavior of (fe near the 
boundary a: = ±l can be represented by 

( 1 . 20 ) (ft.eiV < MV < for x€{yf,l), 

( 1 . 21 ) for xe{-l,yf), 

where —l<yf<yf <1 satisfy limu(?/^) = 0 , and 


( 1 . 22 ) 

M,eiV=C+M' 


'(7+ 

^(l-a:)-flogU+, 

1 

(1.23) 



^(l-ba:)-flog 

1 


Here A^, and Df^, i = l,2, are constants depending on e such that A^^l, 


Bf 




± . 
i,e 


■ \/ai-\- P 2 and D 


± 

i,e 


y/aieXt g+/32 + Vai+t*2 
y/ aie±*-“+/32 —\/“i+/32 


as e goes to zero. 


In the case of fVi = 1, V 2 = 2, oi = 61 = 1 and &2 = 2, we may solve equation (11.191) pre¬ 
cisely and get the form of ( 11 . 221 ) and (11.231) near x = l and x = — \, respectively. One 
may remark how the values c,t,ai and P 2 affect the asymptotic behavior of near the 
boundary a; = ±l. 

When a^fi (the non-electroneutral case), the asymptotic behavior for the solution 
(pe, rie and Pe of the equation (ira-(fr7i) with the boundary condition dElD is stated as 
follows: 


Theorem 1.5. Assume 0<a<f) and (pQ =(pQ- -^ 6 ^ '/’e £C'°°((—l,l))nC'^([—1,1]) be 
the solution of the equation with the boundary condition and Pe > 0. 

Then 

(i) When 0<e<l and 0 <«:<!, there exists a positive constant depending 

on e and k such that limej^oAe)^) = 0 and 

(1.24) ^ — Ae(«:) <ne(x) <pe(x) < ^ -I-Ae(K), for xS [—l-l-e'^,1 —e'^j. 
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Moreover, we have 

(1.25) 

(1.26) lim sup 


limne(±l) = 0 and lime^Pe(±l) = 
eio eiO 


(a-/3)2 


sup 

f \ ^ 

= lim 

sup 

PeW - 0 


z 

£ 4,0 3,g| 


2 


- .= 0 ’ 

p — pi 

(1.27) lim / neix)dx = lim / nJx)dx = Q, 

4oy_i 

— ^ Ol 

Pe{x)dx = lim Pf^{x)dx =—-—. 

1 2 

(ii) Let K be any compact subset of (—1,1). When 0 < e ^ 1 is sufficiently small, the 
asymptotic expansion o/^^(a::) — 0e(±l) in e with the exact leading-order term 
log^ and second-order term 0(1) is given as follows: 

(1.29) (/)e(x)-(()e(±l) =logL-|-log +0£(1), for xGK, 

4a 

where 0 ^( 1 ) denotes as a small quantity tending to zero as e goes to zero. 
Similar results also hold for 0< ft < a. 


Remark 1.3. 

(i) To exclude the boundary layer of with thickness (cf. Theo¬ 

rem 1.6 of m), we consider integrals of n,, and Pe over the in¬ 
terval [—l + e'‘,l~e^]) where 0<«;<1 is independent of e. Note 

that Tie and can be represented by (see and that The¬ 
orem \1.5\f ii) implies rie(x)(ia; == a > 0, 

lime4,o ( n^{x)dx = 0 and lim^o ( Pe{x)dx = 
/3 —a>0. This shows that as e approaches zero, both the total concentrations of 
anion and cation species in the bulk [—l + e'‘,l~£^] tend to the same positive 
constant a, while the total concentrations of anion and cation species in the 
region [—1, —l + e'^)U(l —e”,l] ( which is next to the boundary with thickness 
2e^) tend to zero and positive constant (3 —a, respectively. 

(ii) We want to emphasize that Theorem \1.5Y ii ) improves the asymptotic behavior 
o/(^e(a;) — (/)e(±l) shown in Theorem 1.5 of our previous paper fSOf . 


Following results play important roles throughout this paper. 

(a) Multiplying the equation (O by > (HU may be transformed into 


Ni 




OLk 


(1.30) 


felt j\e°^’^<l>dy)dy 

+ Oj, 


( 31 ) 


Pi 


N 2 

_ 

1^1 J^j^e-^Ne(y)dy 


„-bi4>^{x) 


where is a constant depending on e. 

(b) Differentiating (11.21) to x and multiplying it by (j)(, 

{x)(j)({x) 

(1.31) 

/ 2 ^2 ,20 

^ y- QfcO^fc I ^-bi4>Ax) 
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The rest of this paper is organized as follows: The proof of Theorems 11.31 and 
11.41 are shown in Section [2] In Section [3l we compare the CCPB equation (11.211 and 
the PB equation (11.41) , and give the proof of Theorem 11.11 In Section 31 we consider 
the non-electroneutral case and give the proof of Theorem 11.51 In Section [5l several 
numerical experiments results of the CCPB equation (EH and the PB equation ca 
are presented. The numerical computations are basically preformed using finite element 
discretizations. In the final section, we state the conclusion. 

2. Electroneutral cases: Proof of Theorems 11.31 and 11.41 Let be the 

solution of the equation (11.21) with the boundary condition (11.51) . A crucial property of 
(pe is given as follows: 


Proposition 2.1. Let G l,l))nC'^([—1,1]) he the solution of the equation 

S1.2\) with the boundary condition 11.51) . Then the following properties hold. 

(i) Either 0' has at most one zero in [—1,1], or 0'=0 on [—1,1]. 

(ii) If is nontrivial (i.e., nonzero solution), then 

(2.1) 4>'f{x2)4>'f_ix2)>4>'fixi)(j)',,{xi) for -l<Xi<X 2 <l. 


Proof. We prove (i) by contradiction. Suppose there exist j/i, ?/2 G [~1; 1] such that 
2/1 < 1/2 and (yi) = ())((1/2) = 0. Then integrating (|1.31l) from yi to 2/2 and using inte¬ 
gration by parts, we may get 




•y\ 


/ Ni 

E 


alat 




^ak<t>e{x) 


N 2 

-E 


blPl 




^-hMx) (j)'^^(^x)dx = 0, 


which implies (/)( = cj)'f = 0 on [yi , 1/2]. Here we have used the hypothesis (i/i) = ^!>((2/2) = 
0 and each ak,(3i>0. On the other hand, the CCPB equation (11.21) has the following 
form 


Ni 


N 2 

/=i 


for xG (—1,1), 


where A^^e’s and Bi^fs are constants, therefore pe. satisfies the unique continuation 
property. Therefore, <()( has to be identically zero on [—1,1]. This completes the proof 
of Proposition 12.ll iL 

To prove (ii), we assume that is a nonzero solution of (1121). Thus, for any 
subinterval (a:i,X2) C (—1,1), Proposition 12.ll il immediately implies 


' Ni 


aluk 


W ] (j)f(x)dx > 0. 


Integrating (11.311) over the interval {xi,X 2 ) and using (12.21) . we obtain (12.11) and complete 
the proof of Proposition 12.II □ 


The following interior estimate of (pe is a key step for the proof of Theorem 11.31 


Theorem 2.1. Under the same hypotheses of Theorem \1.3\. we have 

(i) —(/)e(—1) = (/)e(l) > 0 and ^((1) = </>((-1). The solution pe is monotone 

increasing on [—1,1], concave on (—l,x*) and convex on (a:*,l), where x* G 
(—1,1). Moreover, we have 

(2.3) 


rnax I </>, (x) ]< (/)e(1) < </>]]■. 

xG|—1,1] 
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(a) There are positive constants Ci and Mi independent of e such that for 
x£[— 1 , 1 ] anrf 0 <e<Cl, 


(2.4) 



Mi(l-x) 


)■ 


Remark 2.2. 

(i) Replacing ife by </'£ + C' for any constant C, the equation 11.^1) is invariant. 
Hence Theorem \2.1\ (i) implies that for any (j)Q and (pQ , (pe is monotonic on 
[- 1 , 1 ]- 

(ii) When A^i = A^ 2 ; cxk=^k cind ak = bk for as for Theorem 1.2 

in the solution of and mM) is an odd function on [— 1 , 1 ], and all 

denominator terms of (m become equal. Then one may follow the argument 
ofJWf to get the asymptotic behavior of ’s. However, as Ni ^N 2 , otu ^ Pk or 
OkT^bk for some k, the solution p^^ may not be odd on [— 1 , 1 ] so the argument 
of ]3(Jl may fail for this case and we have to develop a new argument to 
prove Theorem \2.1l 


2.1. Proof of Theorem 12.11 

Integrating (fL^ over (-1,1) gives J^iP'f{x)dx = J2k=i^kak-Ylf^ibiPi = 0- This 
implies p'^{l) = p'^[—l) and there exists x*s(—1,1) such that ((>"(x*) = 0. Along with 
the boundary condition (11.51) . we find —^e(l) = ^e(—1). Setting xi=x* and X 2 = x* in 
respectively, we get 

(2.5) 0"(x)(}()[,(x) > 0 for X € (x*, 1], and ^"(x)(/)((x) < 0 for x G [—l,x*), 

which implies: (a) Both p'^ and p'f never change sign and share the same sign on (x*,l]; 
(b) p'^ and p'f never change sign and have opposite signs on [—l,x*). Consequently, (/)( 
never changes sign on [—1,1] due to (a), (b) and ()i((l) = ()i((—1). Now we claim <()( > 0 
on [—1,1]. We state the proof using contradiction. Suppose p'^<Q on [—1,1], then the 
boundary condition (11.51) implies = 0e(—1) — ? 7 e(/)((—1) > 0e(l)+ 7 ye^((l) = , which 

gives a contradiction. 

Therefore, we get </>(:>0 on [—1,1]. Along with the boundary condition (II.5p . we 
prove (lOl) . 

Furthermore, by (12.51) and p'^{x) > 0, we have p'f (x) < 0 for x G (—l,x*) and p'f (x) > 0 
for xG (x*,l). Hence we complete the proof of Theorem 12.II fib 
By (12.31) and (11.311) . we obtain 

(2.6) e^{<P?{x))" > 2e^p'f'{x)P',{x) > AMIP'^{x) 

1/2 

for xG (- 1 , 1 ) and e> 0 , where Mi = \ alakc'^^'^^o bfPie^^^'f’o^ . 

Note that (/)((1) = </)((-1) >0. By (|2.6p and the standard comparison theorem, we get 

(2.7) 0<p^{x)<p^{l)(e - +e ^ j. 

It remains to deal with ^^((1). By (12.31) . there exists x^G (—1,1) such that 0<(()((xe) = 
^ 0 +^ Subtracting (|1.30l) at x = Xe from that at x = I and using (12.31) . it is 
easy to get </'((!) <^ as 0 <e<Cl, where Ci is a positive constant independent of e. 
Along with (12.7|) . we get (12.41) and prove Theorem 12.H jib 
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Therefore, we coniplete the proof of Theorem 12.II 

Note that (|1.30l) plays a crucial role on the asymptotic behavior of as e 4,0. The 
estimate of the constant in (11.301) is given as follows: 

Lemma 2.3. Under the same hypotheses of Theorem \2. 11 we have 

(i) For any x,y€ (—1,1), (j)^(x) — 4>^{y) converges exponentially to zero as e goes 
to zero. 

WJVi JV2 \ 

(a) limCe = —- [ ) 7 where is the constant defined in hl.30\) . 

^ \fc=i ) 


Proof. (12.41) implies that for any x,y G (—1,1), 
lim^^(a;) = 0 and 


( 2 . 8 ) 


etO 


|A(a:)-A(2/)l<;^( 


e ^ +e ^ 


Ml (l + iy) Mi{l — y) 

e ' +e 


This may complete the proof of Lemma lOT i). Note that (1^4^ gives 
sup„, i)|A(a::)-A(2/)l<4C'i/Mi and lim^o I - A(y)l = 0 for x, j/e (-1,1). 
Applying Lebesgue’s dominated convergence theorem, we obtain 


(2.9) 


lim ■ 


afcC 


aicipUx) 


^10 J dy 2 


= and lim ■ 


Ae 


-bi4>c{x) 


A 


'^ 4 '^ J_ie~^Ue(v)dy 


for k=l,---,Ni and Z = 1,• • •,V 2 . 

Therefore, by (ll.30p . (12.81) and (12.91) . we prove Lemma HSlii) and complete the 
proof of Lemma [2.31 □ 

2.2. Proof of Theorem 11.31 To prove Theorem 11.31 we need the following 
lemma: 


Lemma 2.4. (i) Under the same hypotheses of Theorem \2. 1[ we have 


{ii) If rje ^ 0, then 


Ni 


( 2 . 11 ) 


oak4>t{i) p^^-bi4i,{i) 




2t), 


^/_^e“'=‘Z'Ay)dy ^ J_j^e-^‘l’-(y'idy 


Proof To get (12.101) . we subtract the equation (11.301) at x = — l from that at x = l. 
Here we have used the facts that (^'(1) = (()'(—1) and (j)e{—l) = —(j)ei^) which come from 
Theorem 12.11 (i). Setting a: = l in (11.301) . we use ()1.5|) to get (I2.11|) . and complete the 
proof of Lemma [HH □ 

To uniquely determine the values c and t, we need the following lemma. 

Lemma 2.5. Assume'Y^^f.^akOLk = 'Yl!i=i^iPi- Then 

(i) f is strictly increasing on (0,oo) and strictly decreasing on (—oo,0). 
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(ii) There exists a unique solution (t,c) of the equations S1.16\) - n.lS\) . 

Proof. By (11.111) and akOLu = hPi, it is easy to check that f'{s) >0 if 

s > 0 and f'{s) < 0 if s < 0. This shows (i). To prove (ii), we need 
Claim 1. There exists 0 < s < —27(/(s) —/(O))^/^ such that 

f{s) = f(^s-2<j)^ + 27(/(s)-/(0))^/^) . 

Proof. [Proof of Claim 1] Let fc(s) = s —2^0"+ 27(/(s) —/(O))^/^ for sSR. Then 
k{0) = —2(l)^<0, k{2(l)^)>0 and fc'(s) = 1+ 7(/(s) —/(0))“^/^/'(s) > 0 for sG(0,oo). 
Hence, there exists si G (0,2 ^q ) such that fc(si) = 0 and 

(2.12) k(s)<0 for sG(0,si). 

Let h{s) = f{.s) — f{k{.s)) for sGM. Then /i(0) =/(O) —/(—2 (()q ) <0 and h{si) = f{si) — 
/(O) >0. Hence there exists S 2 & (0,si) such that h{s 2 ) = 0. On the other hand, (I2.12p 
implies fc(s2)<0i.e. S 2 <2(j)Q - 2 j{f{s 2 ) - fiO))^/"^. 

Therefore, we complete the proof of Claim 1. □ 

Now we want to prove Lemma l2.5l iiL By Claim 1, 

(t,c) = - 7(/(s) - /(0))l/^(/);j■ - s -7(/(s) - /(o))^/^) 

is a solution of (jl.l6|) and (jl.l7l) . Moreover, 0 < s < 2(j)Q — 2j{f{s) — /(O))^/^ gives jcj < 
t<(l)Q. Hence (I1.16|) - (I1.18|) have a solution. The uniqueness of (ii) can be proved by 
contradiction. Suppose (ti,ci) and (0,02) solve (II. 1611 - (II. 181) and ti>t 2 . If f{ti — 
ci)> f(t 2 -C 2 ), then -t 2 > (/>[[ - ti =7(/(ti - ci) - /(O))^/^ > 7( 7(0 - C 2) - 7 (0))^^^ 
i.e. </'[)’ —12 >7(/(t2 — C2) — /(0))^/^ contradicts (t2,C2) a solution of (I1.16I) - (I1.18I) . Thus 
/(ti — Cl) </(t 2 — C2) and then (11.171) gives /(—— Ci) </(—^2 — C2). Furthermore, by 
Lemma I2.5f ii , we obtain ti — ci < t2 ~ C2 and —fi — Ci > —^2 — C2 which implies 11 < ^2 a 
contradiction to the hypothesis H >^2- Hence, ti =t 2 :=t*. Here we have used the facts 
that 11 — Cl, ^2 — C2 > 0 and —11 — ci, —^2 — C2 < 0. 

To prove ci=C 2, we set g{s) := f{t* — .s) — f{—t* — .s). Note that Lemma [13] (i) 
implies g'{s) = — f'{t* — s) + f'{—t* — s) < 0 for jsj < t*, i.e., g{s) is strictly decreasing on 
{—t*,t*). Therefore, we have ci =C2 and complete the proof of Lemma [231 (ii). □ 

Now we shall give the proof of Theorem 11.31 


Proof of Theorem 11.31 By Lemma [2.3f ii. it suffices to prove lim(/)e(0) = c. By 
Theorem 12.11 {|(/'e(0)|}^^Q has an upper bound. Then we set limsup(/>e(0) = Cs and 

£4,0 

liminf(/)e(0) = Ci. Hence there exist sequences {cilieN and {cilieN tending to zero such 
that lim (j)^ {0) = Cs and lim ())j- (0) = Ci. We may rewrite (12.101) and (I2.11|) as follows: 

j—¥oo 


Ni 

E 


Oik (1 


3afc(0e(l)-0e(O)) _ p-afc((^e(l)+<;i. 


(0))^ 


gak(<t>e(y)-<f>e(0))^y 

Pi _ g-bl{4>sp)-<l>e{0))^ 


N 2 


=E 


1^1 


(2.13) 


































Boundary layer solutions of Charge Conserving Poisson-Boltzmann equations 


14 

and 


(2.14) 


g2 ^ afc(0,(l)-0,(O)) 

^ O -bi(,pAl)-4>e(0)) 

_ \^(J 


We divide the proof into two cases. 


Case 1. Iime 4 ,o^=oo. 

Note that |())g — ())e(l)| < 2 (()q . By (12.9L (ll.ll|l . (12.1411 and Lemma 12.31 (ii), 
we have /(limsupe|o(())e(l) - <(>e(0))) = /(liniinfe 4 ,o(())e(l) - i^e(O))) = /(O) - Then 
by Lemma 12.51 (i), lim(: 4 ^o(</>e(l) — ^e(O)) = 0. Along with (12.131) . we find 
/(-hm 4 o((('e(l) + </'e(0))) = /(0), this gives lime^o(«f>£(l)+ </'e(0)) = 0. Consequently, we 
have lim(: 4 ^o<(>e(l) = lime4,o<(>e(0) = 0. Hence, we obtain c = t = 0 and complete the proof 
of Theorem 11.31 lib 


Case 2. limcj^o^ = 7 < oo. 

By Theorem 12.11 has an upper bound. Then there is a constant ts and 

a subsequence of {e^} (for notation convenience, we still denote it by {cj}) such that 
hmj^.oo0e^(l) = ts. Putting e = €j in (12.131) and (12.141) and using Lemma [2?3r iiL one 
may check that {ts,Cs) satisfies 

(2.15) {(j)^-tsf=-/^{fits-Cs)-f{0)) and f{ts-Cs) = f{-ts-Cs). 


Now we claim 


t/gA^ 0, then we have |c. 


S I ^ ^ 


|Cs| Since |(/>£(0)| <(/>(!) < (/>d 

9 ^ 0 . If |cs|=ts, then by the second equation of (12.151) and Lemma [23r il . we have 
ts — Cs = —ts — Cs = 0, i.e. ts = Cs = 0. Along with the first equation of (|2.15l) we find 

\<ts<4>0- Along with (|2.15p . its,Cs) 


4>q =0, which contradicts to ^g 7^0- Hence \cs 
satisfies (|1.16ll - (|1.18p . 

Similarly, there is a positive constant ti such that {ti,Ci) satisfies (I1.16ll - (ll.l8p . By 
Lemma [2.51 liil. we get Cs = Ci = c and ts=ti = t, where limej^o(0) = c, limej,o</'e(l) 
and (t,c) satisfies (I1.16|l - (ll.l8p . Therefore, we may complete the proof of Theorem 
0(11). 

By (II.up and |t — c| < 2 (()q , f{t — c) — /(O) is uniformly bounded for all 7 > 0. Conse¬ 
quently, by (I1.16I) - (I1.18L wehavehmg,_>ot(7) = (()g and hmg,_>o 0(7) = c*, where |c*| <4>q 
is uniquely determined by /((/ig — c*) = f{—4>Q — c*)- By (|1.16l) we have f(t — c) — /(O) = 

) which and (11.171) give limg,^oo/(i — c) =lim..y^oo/(—i — c) =/(O). By Lemma 
I2.5l i') and the continuity of /, we find limg,^oo(^ —c) = limg,^oo(—t —c) = 0. Hence, 
lim.y_>oo^ = hmg.^ooC = 0 and complete the proof of Theorem 1 1. 3l ii- A1. 

It remains to prove Theorem II.3l ii-Bl. By (II.16p - (II.18p . t and c are uniquely de¬ 
termined by 7. Hence we can consider t = t{'y) and 0 = 0(7) functions of 7. Due 
to /(t(7) —0(7)) —/(O) 7^0 on (0,oo) (by (11.181) and Lemma [2^ iH. t{'j),c{-y): (0,oo)—)> 
{—(j)Q,(j)Q) are continuously differentiable. Differentiating (11.161) and (|1.17l) to 7, one 
may check that 


(2.16) 


l- 7 (/(i-c)-/( 0 )) 


nt-c) 


f'{t-c)-f'{-t-c) 


dt 

c?7 


= -(/(t-c)-/(0))i/^ 
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and 

(2.17) 


d 2f'{—t — c) dt 

d'y^ ^ f'{t-c)-f'{-t-c) dj' 


If ^ changes the sign on (0,oo), then there is a 7*g(0,oo) such that ^( 7 *) = 0. By 
(|2.16l) and Lemma [2?5l' i'). we have ^( 7 *) = 0 ( 7 *) contradicting to (|1.18|) . Hence ^ keeps 
the same sign on (0,oo), and then Theorem ll.3f iiii gives ^ < 0 on (0,oo). On the other 
hand, Lemma IS^T i') and (11.181) imply that both — and — are positive. 

Consequently, by (I2.17|) . we obtain that -^{t — c) and ^ share the same sign. Therefore, 
we prove Theorem 11.31 (ii-B) and complete the proof of Theorem 11.31 


Remark 2.6. Suppose limgj^o —= 0. Then Theorem \2.1\ (i) and Theorem \l.S\ (i) give 
>-0 uniformly m [—1,1] as e],0. 

2.3. Proof of Theorem 11.41 

For convenience, setting v{x) = (j)e{x) — c, By (11.301) . we find 


2 

(2.18) —v'^{x)=Y^ 

k- 


Olk 


„akv(x) _ 


N 2 


Pi 




f^J^_^e-bMv)dy 


^-biv{x) 


+ Cf. 


Note that by Theorem I2.11 ii and Theorem 11.31 we have |r’(a;)| <</>({’ +|c| and 
lim / = lim / e~^‘'^^^^dy = 2. Hence by Lemma [2.3r iiL it is easy to check 

that 

(2.19) \e^v'‘^(x) - [f{v{x)) - /(O)] | < (5(e), 

for all xG [—1,1], where (5(e) is a positive quantity tending to zero as e goes to zero. By 
Theorem 12.If il and Theorem ll.3l iiL we have v’ = (j)'^>0 on (—1,1) and 


(2.20) limti(l) = t — c>0> —t — c = limu(—1). 

e^-O €4,0 

Thus there exist e* > 0 and ?/+€(—1,1) such that v{yf ) = 

— log|l— [5(e)/(ai +,82)^] > 0 for 0 < e < e*. Note that limej^oi')?/^) =0 and 

(2.21) u(a;)>-y(j/+) = -log|l-[(5(e)/(ai+/32)^]^^'‘|>0, Va:G(y+,l). 

Now we begin to deal with (12.191) when iVi = 1, N 2 = 2, ai — bi = l and &2 = 2. Here 
we have Oi = Pi+ 2(32 and f{v{x)) — /(O) = (1 — +P 2 )- Note that such a 

formula is valid only when 7Vi = l, N 2 = 2, ai=bi = l and 62 = 2. Along with (I2.21|) . we 
have 


iG(e)(l-e-''("))2(aie''(")+^2)+<5(e)</(u(x))-/(0)<(A-e-"("))2(aie"(")+;92)-<5(e), 

( 2 . 22 ) 

for xG (j/^,1) and 0 < e < e*, where K{e) = l — \/S(e) and = 1 +i/S(e)/(ai + P 2 )- Con¬ 
sequently, (|2.19I) ~ (I2.22I) give 

v'(x) ^ 1 ^ 1 v'(x) 

(t^-e-^(=^l)x/aie’'(^^+P 2 ~ ^ ~ \/K{e) (1 - e-"(^))^Qfie“(®)-b /32 ’ 


(2.23) 
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for xe{y+,l) _ 

Integrate (12.2311 over (y,l) for yG (y+,1), we obtain 


v\x) 


(2.24) 


(„_g-«(a::))v/Q-ie«(^)+/32 


dx 


1 {\/ +/32)a- \/ai +/32a)(\/+/32)a +Vai +/32C 


y aia + /32a2 (\/ (Q!ie"(i) +/32)a + \/ai +/32a)(\/(aie"(^) +/32)a — -v/oi +/32a) 

for a>0. Hence, (12.231) and (12.241) imply 

<e’'W 
C’*’ 

— (l-a^)+log U+ 



) csch^ 

'(j+ 

—^(l-2:)+logU)*;^ 

V aiy 

/ 



(2.25) 



<i+( 

for 

a; e (2/^,1), 

where 


y/ aie” 

('-W/32 + \/“lte+)92 

and 

n+ —: 

^2,e — ^ 

V aie” 

('-W/32 —\/“lte+/32 


cscm 


C+=Kie)V^T+^, D+ = 


lime;oie=lime.to-?^(e) = l, we get (jl.2b|) . 

Similarly, we also have (jl.21l) . 

Therefore, we complete the proof of Theorem ll.4l 

When 7 Vi = V 2 = 2, ai = bi = i, i = l,2, and oi+ 2a2 =/3i + 2/32, we may follow the 
similar proof of Theorem 11.41 and obtain the following result. 


Corollary 2.7. Under the same hypotheses of Theorem \2.1l suppose Ni = N 2 = 2, 
ai = bi = i, 1 = 1,2, and ai +2a2 = /3i + 2/32. Then 


(2.26) 

(2.27) 

where 


<4>e(.x)<(t)'^^^(x), ^XG{Xe,l), 
< </>e(a;) < <(’^e(a:), Vx G (-l,Xe) , 


4 'te{x)=C + \og 


coshhfg(x)± 




cf 


Keix) = -^(1 =F x) +logi7i „ / = 1,2. 


coshh±(x) + ^ 

Here A=l + ^, 5=^(1 + ^) and Cf^/s, H^/s, / = 1,2, are positive con¬ 


stants depending on e such that 


limC±=ya2[(4l+l)2-H2], 

e4-0 ’ 


ef,0 ’ 


± t —c 


A-B-\-e 
A + B + e^^-’^ 


A-B-GI 

a+b-gi 


^ jA-B + e^*-^^ lA-B + 1 


A + B + e 


i t—c 


A + B + l 


3. Proof of Theorem 11.11 In this section, we study the asymptotic behavior 

of solution (/)e of the PB equation (11.41) with the boundary condition (11.51) and give the 
proof of Theorem ll.il Surely, the PB equation (11.41) can be transformed into 

e^(l)'f{x) = ^f{Mx)), 


(3.1) 
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where f{s) + is defined by (11.111) . It is well-known that the 

equation (11.41) has the unique solution (j)^ £ C'°°((—1,1)) (~) C^([—1,1]). As for (11.301) . we 
use (EU to derive the following identity 

(3.2) = + 

Moreover, we use the similar argument of (I1.31I) - (I2.1I) to get 

(3.3) (j)"{x 2 )(l>e{x 2 )>(l>”{xi)(j)'^{xi) for—1 <xi <a;2 < 1. 

Using standard maximum principle to dUD and CH), we obtain 

(3.4) min{0,(/)J,(;iQ } < (/)e(x) < max{0,(;i;[,())o }, 
for xG [—1,1]. 

Now we state the proof of Theorem 11.11 

Proof of Theorem 11.11 Multiplying the equation (11.41) by (f>e, we obtain 

/ Ni N2 \ 

(3e§(j(/)^(a:))">2e^(/)"(a;)(;ie(a:)= I | </>£(x) > Cs(x) 

V/c=l 1=1 ) 


/ JVi N2 \ 

where C^= inf s“^ I Yafca/ce“'“® — I >0. Here we have used the hy- 

sGR.s^O \ ^^' / 

\fc=i 1=1 / 

pothesis Y^^Li^kOik = Y^f=ihPi to assure as a positive constant. Thus by (13.41) . (13.Sp 
and the standard comparison theorem, we get 

(3.6) |(/)e(x)| <max{|(/)+|,|(/)(C|} , VxG(-1,1), 

and complete the proof of Theorem ll.ll il . 

Suppose min{())Q ,(^Q } > 0. Then (13. 4p gives 0 < (()e(x) < max{(/)(}",<))((■}, together 
with dSD) and Lemma E5](i), we may find </>[(> 0 on [—1,1]. Here we have used the 
hypothesis that Y^kLi^kO-k = Y^^=i^iPi- To complete the proof of Theorem ll.il fii. we 
need to claim: 

Claim 2. There exist e*>0 and x*g(— 1,1) such that (j)Ax*)= min (jiAx) for 

a:e[-l.l] 

0<e<e*. 

Proof. We state the proof of Claim 2 by contradiction. Suppose </)( preserves the 
same sign on (—1,1), Ve>0. Without loss of generality, we may assume (f'^ix) >0 for 
xG (—1,1). Then by (II.5L one may get 4>e{x) > (()e(—l) ></'(]' > niin{(/)(]■,(()((■}. Along with 
(13.6p . we obtain O = lim£j^o0e(O) Y}> '"^hich contradicts to the assumption 

min{^Q ,0Q } > 0. Consequently, there exist e* > 0 and x* G (—1,1) such that (/)((x*) = 0 
as 0<e<e*. As for the proof of Theorem 12.11 lii. we may use (13.3p and the fact that 
4 >' f >0 on [—1,1] to get (()((xi) < 0 < ((i((x2) for xiG(— l,x*) and X2G(x*,1). Hence, 
(pe attains the minimum value at an interior point x*g(— 1,1). This completes the 
proof of Claim 2. □ By Claim 2, we complete the proof of Theorem ll.il (ii). Similarly, 
Theorem o (ill) can also be proved. 

We prove Theorem 11.11 (iv) in two cases: (I) (p'l never changes sign on [—1,1]; 
(II) (p'f changes sign on [—1,1]. For the case (I), without loss of generality, we may 
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assume >0 on [—1,1]. Then by m, 4>e >0 on [— 1 , 1 ] and the maximum value of 
(j)e occurs at the boundary x = ±l. Suppose (/)e(l) = maxa,g[_i_i]0e(x). Then (/>'(!) >0. 
Moreover, by the boundary condition m, we get (j)Q = > 0, which gives 

(/)q <0 due to min{(/)Q ,(/)q } < 0. Consequently, = — Hence by 

the assumption of on [— 1 , 1 ], we have 4>'^{x)><j)'^{—l)>0 for xG[— 1 , 1 ], i.e., 

is monotone increasing on [—1,1]. Along with the boundary condition (11.511 . we have 
<^0 < < 4'e{x) < < (l>o for X £ [-1,1]. Similarly, if (/)e(-l) =max„,g[_i_i] (/)e(x), 

then we obtain (/ig < 0 on [—1,1] and (j)Q < ipeix) < 4>o ■ This proves (hi). For the case 
(II) there exists Xg G (—1,1) such that (xg) = 0. Then we may use (13.31) and the same 
argument as in Theorem 12.11 (i) to get Theorem 11.11 (iv). 

It remains to prove Theorem EHv). Using (IT^ . Theorem ll.ll il and the similar 
argument of Lemma 12.31 we may obtain hmgj,oC'( = — ^/(O) which implies 

(3.7) nm[((;iJ-(/-g(l))2-72(/((^g(l))-/(0))] =0, 

by setting x = 1 in (]3.2I) and using boundary condition (II. 5p with limgj,o — = 7 . There¬ 
fore, as for the proof of Theorem 11.31 livl. we may use Theorem [O] (i)-(iii) and (13.7|) 
to get Theorem ll.H vl and complete the proof of Theorem 11.11 

Remark 3.1. If have limg 4 .o'iu(a;) = r for all xG(— 1,1), 

where r is uniquely determined by f'{r) = 0. The proof is similar to the proof of 
Theorem 4-2 ofm- 


4 . Non-electroneutral cases: Proof of Theorem 11.51 

In this section, we assume 0 < a < /3 and (j)Q =(j)Q . To prove Theorem [LSI we 
need the following properties, which can be obtained from [30) . 

(PI) Gradient estimates of (()g (cf. Theorem 3.1, [50]): The unique solution i^g is 
even and satisfies 4>”<0 on [— 1 , 1 ], and 4 >'^{xi)> 0 ><j)'^{x 2 ) for xi£[— 1 , 0 ) and 
X 2 G( 0 , 1 ]. Moreover, satisfies 


(4.1) 
and 

(4.2) 


-0'g(-l) = 0((l) 


a — (3 


< 0 , 


m^)\< 


13 —a ( v^(i+x) 


+ e 


V 2 ^(l-x) 


). 


VxG(-1,1). 


(P2) Interior asymptotic behavior of i^g (cf. Theorem 1.5, [5^): For any com¬ 
pact subset K of (—1,1), there holds 


(4.3) 


sup 

0<e<l 


(()g(x)-(()g(±l)-log 


1 


< 00 , VxGAT. 


(P3) Estimates of Ug and p^-. In [30], we have established the following estimates 
(see (3.9), (3.15) and (3.37) of [50]1: 


(4.4) 


A<J e^^^y^yj e-‘^‘^yUy<^, 


ae 


ddo) 


- 0 .( 0 ) 


^2 /-i 




f>?{y)dy = 


a + f3 


(4.5) 


-1 


2 











C.-C. Lee, H. Lee, Y. Hyon. T.C. Lin and C. Liu 


19 


and 

(4.6) 


{a-l3f 

8e2 


< 


5e ^ {a — PY , ct + P 


J^^e-‘l>-^y')dy 8e2 


Using (11.71) and the fact that p^{l) = p^{—l) and p'^{0) = 0 (by (PI)), we can 
transform (ioi)-(ra into 


(4.7) 

(4.8) 

(4.9) 


O" , /X / X , C(P 


^ <neix)peix)<—, Va;G[-l,l], 


n<:(0)+p<,(0)- 


^2 ri 


a + P 


(«-/3)^ 

8e2 


<Pe(l)=Pe(-l)< 


' iy)dy = - 2 > 

{a — PY a + P 


8e2 


respectively. 

Having (P1)-(P3) at hand, we are now in a position to prove Theorem 11.51 


Proof of Theorem 11.51 Let /ef= = [—l + e'^,1 —e'^], where 0<e,K<l. For any 
y G le-x, we may use (14.21) to get 




> — a 


n 


\/2a(l + x) \/^(l-x) ' 

e 2= -pe 2. 


dx 


^2V2(P-a) 

<- — -e ^ 

y/ae 

As a consequence, we have 

(4.11) \p^{x) - pe{y)\< \(t)e{x) - ^!)e(0)| + \pe{0) “ (t>e{y)\< ^ 

y/ae 

for x,yGle^. Note that hmej.o e~=0 for 0 <k< 1. Thus (14.111) gives 


(4.12) 


lim sup \peix)-peiy)\=0- 

X,yGle>^ 


For 0 < e < 1, we may set x = 0 in f|4.3p and combine the result with f|4.1Qp to get 


(4.13) 


sup 

0<e<l 


p{y) - p{±l)-log 


1 


<00, WyGle'^. 


To prove (ll.24p - (ll.26p . we need the following claim: 


Claim 3. 

(i) At the boundary a: = ±l, we have 

(4.14) lim^^^^ = 0 and limeV(±l) = 

^ ^ 40 £2-^ eto ^ 8 ’ 

for any r > 0. 

(ii) Assume 0 < e < 1. Then there exists Ae(K) > 0 such that limcj^o K{n) = 0 and 

/,x_x a , , X /X y/ aP a , ^ P ^ ^ 

(4-15) -:^-Xe{K)<ne{x)<^^, <Pe(a:)< | + A,(k), 

for ccG lei^. Moreover, 

(4.16) lim sup |ne(x) —ne(0)| =lim sup |pe(a:) — pe(0)| = 0. 

eio ffo xei^K 
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Proof. dUZl) and (gj]) give < ^^(-1) =ne(l) < This shows 

lime^o = 0 for any r > 0. Along with (j4.9p . we prove (14.141) . 

By (PI) and (II.6|) . we have 

(4.17) (()e(0)= max (l)e{x) 
and 

(4.18) ne{x)—pe{x) = e^(l)”{x)<0, Va;S[—1,1]. 

Along with (14.7p . we obtain 

(4.19) Pe{x)>^ and n^{x) < , Va;G[—1,1]. 

By (11.71) . (14.81) and (|4.17l) . one may check that 

(4.20) 0< n,(0) - Ueix) = n,(0) (l -(^1 -e-/^d^)-0do) j ^ 
and 

(4.21) 0 < p,{x)-peiO) =Pe{0) (^e-'^dx)+M0) _ j'g-0.(x)+0.(o) _ ^ 

Consequently, by (I4.12L (14.201) and (I4.21|l . we get (14.161) . 

It remains to prove (14.151) . Let 


(4.22) Ae(K)=max'| sup \n^{x) — nc{0)\, sup |pe(a;) —pe(0)| ]• >0. 

X^I^K 

By (I4.16p . we have limg4_oAe(/c) = 0. Using (14.171) . one may find 


(4.23) 

and 

(4.24) 


= (0) = 7T 


> 






/3e 


-0e(o) 


/3 


J_^e-‘i>dv)dy J_^e-'^dy)+Mo)dy' 


< 


Hence, (I4.19|) and (14.231) immediately give > ne(a:) > n(:(0) — Ae(K) > ^ — Ae(K), for 
xG On the other hand, by (14.191) and (14.241) we obtain ^ <Pe{x) <Pe(0) + Ae(K) < 
^ + Ae(K), for xGI^k. Therefore, we get (I4.15|) and complete the proof of Claim 3. □ 
(11.241) immediately follows from (|4.18|) and (14.151) . and p.25p follows from (14.141) . 
To prove (|1.26p, we rewrite n,{0) = 


(4.25) n,(0) 

By (|4.17|) . we have 



a 

+ e‘^dy)-4>do)dy 




(4.26) 


0< 


eMy)-Mo)dy<2e'^. 
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On the other hand, by (14.121) we get 


(4.27) 


/■l-e" 

lim / = 2. 

elO 


Combining ()4.25|) - (I4.27|) . we conclude that 


(4.28) 


a 


limrif (0) = —. 
€ 4.0 2 


To deal with the limit of value ^£(0) as e tends to zero, we need the following 
estimate: 


(4.29) 


P€(0)-f 


<|ne(0)-p,(0)|+ n,(0)-- 


a 


<\neix)-peix)\+ n^iO)-- +2 X^{k), \fxGle‘ 


Here we have used (14.161) and (14.221) to get the second line of (I4.29|) . On the other hand, 
by integrating dm over /jK and using dm and (14.211 . we obtain 


(4.30p< 




{Peix) — ne{x))dx = + e”) - ((>((1 - e'")) < 4(/3- a)e . 


/-l+e» 


Note that 0 < /c < 1. As a consequence, by (14.2811 - (14.301) we find 


(4.31) 


limp.( 0 ) = -. 


Then (11.261) follows from (14.161) . (14.281) and (|4.31|) . 

By (|4.19p . we immediately get (11.271) . Now we shall prove (|1.28l) . Let g{x)€ 
C'^([—1,1]). Multiplying (11.611 by g{x) and integrating the result over (—1,1), we have 


{n^{x)-pe{x))g{x)dx = 


(l)”{x)g{x)dx 


(4.32) 


a — P 


{g{-l)+g{l))-e'^ j ^(j)'^{x)g'{x)dx. 


Here we have used the intergration by parts and (|4.ip to get (14.3211 . On the other hand, 
by using dm, one may check that 


(4.33) 


(j)'^{x)g'{x)dx 


<{l3-a) max |g(a:)| [ ( 


<2(/3 —a)\/— max |(4(x)| e. 

a \xe[-i,i] 


V^(l + x) \/^(l-x)' 

e 2€ -pe 2€ 


dx 


By ()1.27p . (14.181) . (I4.19L (14.301) . (I4.32p and (14.331) . we have 


+y ^^P4x)ff(3:)dx- ^-^(g(-l)+g(l)) 

) ne(x)g(x)dx + / (ne(x)—pe(x))g(x)dx + e^ / (j)'^{x)g'{x)dx 

d-l+e'' J-i 

sj aPe'‘-|-4(/l —a)e~2ei-~ -\-2{j3 — a) 


n-H-e'' 

-1 


< max \g{x) \ x 

xG[-l.l] 
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Note that 0 < k < 1. Consequently, 

(4-34) 

In particular, let gS C^([—1,1]) satisfy g{x) = 1 for xG [—1,0], g{x) G [0,1] for xG [0,1/2], 
and g{x)=0 for xG [1/2,1]. Then (14.3411 gives lime4,o pe{x)dx= Similarly, 

we have lim^j^o J^_^^Pe{x)dx= Therefore, we get (11.281) and complete the proof of 

Theorem (THi). 

It remains to prove (11.291) . By (11.711 . we have 
(4.35) _ 0^(1) _ log = log 

Pe{x) 

Note that for any compact subset K of (—1,1), we have K Cl^'^ as 0<e<Cl 
is sufficiently small. Hence, by (11.2611 . (14.1411 and (14.3511 . we conclude that 

lim^o {(t>e (x) - (1) - log ^) = log = log uniformly in K. There¬ 

fore, we get (11.2911 and complete the proof of Theorem 11.51 

5. Numerical experiments 

In this section, we do numerical computations to compare solutions of the CCPB 
and PB equations. All numerical results are obtained using the convex iteration method 
[snuniiziiiH] and the finite element methods with piecewise linear space which is used 
to solve the linearized equations. The computational domain and the mesh size h are 
fixed with is [—1,1], h = 2“^^, respectively, throughout the numerical experiments. The 
values of e are set by e = 2“-^ ,/ = 1,3,5, in order to observe the tendency of the associated 
solutions i/e’s as e goes to zero. 

As for [3D] , the numerical scheme can be extended to the CCPB equation (11.21) with 
the boundary condition (11.51) and it can be presented as follows: 


(5.1) 


m+. 


Ni 

yb gfe^fc 


E 


Z=1 




(5.2) (jirn+l —S^^_|_^-|-(l s)(j)rm 

for m = 1,2,'-', where s is a positive constant satisfying 0 < s < 1 with boundary condi¬ 
tions 


(5.3) (/„+i(-1)-?7,^(„_^i(- 1) = 4 , (/^+i(l)-h77,0y_^i (!) = (/+. 

Let i'm+i- =4’m-\-dm with the correction term 5m which satisfies 

(5.4) (5„(-l)-??,(5(„(-l) = 0, 6mii) + ge5mii) = 0 

so that (fim+i = (t>m + s6m = 4>i + 'sSfci linim-j-oo ji^m | = 0, then the iterative scheme 
converges. 

Define the residual function TZ{(j)m) as 


(5.5) 


Ni 


^ ( I 




^ /_ie“'“'^”‘dx 


E 




dx 


— e 


n 

m ■ 
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Then we obtain 

(5.6) = = 

Integrating 'R,((j)m+i)— we may use (I5.5p and (j5.6p to get 

TZ{4>m+i)dx = {l-s) f TZ{(j)m,)dx. 

1 J-i 

In case of s = 1 in (IE3, numerical scheme may not converge and oscillate during the 
iteration procedure. When 0<s<l, we have empirically observed that the value of 
s should be compatible to Ce^ in order to let the iteration converge. Moreover, the 
value of C is chosen in the interval (0,1) so that the convergence of the scheme can be 
guaranteed. In the iteration procedure, the value I0“® is applied for stopping criterion 
with I |(Im I |oo — 11 (0m+l ^m)/'5||oo- 

For the PB equation (II.4L we replace the denominators of the right hand side of 
the equation (EH) by the value 2. Then as for the scheme of (I5.1l) - (l5.3p . we have a sim¬ 
ilar way to solve the PB equation (11.41) with the boundary condition (11.51) . numerically. 
To compare solutions of the PB and CCPB equations, we firstly set the parameters 
as 7 Vi = 1,.^2 = 2, ai = 6i = l,&2 = 2 and oi = 1.2,/3i =/32 = 0.4 so that the electroneu¬ 
tral condition oioi = &i/3i-1-62/32 holds. The numerical computations also impose the 
boundary data as (j)Q = —(t>Q = 1 and the values t/^’s for the boundary conditions m as 
r/e = 0.5e^ and 0.5e which include the cases of lime^o ^ = 0 and 0.5. The corresponding 
results are presented in Figure EH] and Table EH) consistent with Theorem II.31 and ll.II 

In Figure E31 one may see the difference between the solutions of (IL2|) and dLl 
with the same boundary condition dlHI and the valence Zi = —1 for the anion, * = I, and 
Zj = 1,2 for the cations, j = I,2, respectively. The solution profiles of the PB equation 
(11. 4p are plotted as (red) dash-dotted curves and those of the CCPB equation (11.2p 
are sketched as (blue) solid curves. Here the index numbers, 1,2,3 are associated with 
various values of e’s, and a (black) dotted line is represented as the axes for a reference. 

Table ETT] shows the numerical results of (('e(O) and c for the CCPB and PB equations 
where the value c is defined in Theorem ll.3l can be computed by Newton’s method. One 
can easily see that for the PB equation, the value c is always equal to zero but for the 
CCPB equation, the value c may not be equal to zero. The ratio /3i//32 may affect the 
value c and t. As /3i//32 varies, the numerical values of ((>e(0), <('e(l), c and t are presented 
in Table E2] for the case of ai = 6i = l, 62 = 2 and e = 2“®. Note that the numerical 
values of (j)e{0) and 4>e{^) are quite close to those of cS (c*,0) and t, respectively. This 
is consistent with the results of Theorem 1 1.31 We remark that if t is fixed and Pi/P2 is 
decreasing, then the value c is decreasing. 



Table 5.1. The numerical results o/(j!)e(0) and its limit value c of PB and CCPB equation in 
Figure 15171 



e 

2“^ 

2"^ 

2"^ 

C 

I 

PB 

CCPB 

0.0106 

-0.0459 

0.0000 

-0.0964 

0.0000 

-0.1081 

0 

-0.1126 

II 

PB 

CCPB 

0.0079 

-0.0311 

0.0000 

-0.0442 

0.0000 

-0.0442 

0 

-0.0441 
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X 


II. Comparison of cj)^ of PB and CCPB, t)^ = 0.5e, 7 = 0.5 



X 


Fig. 5.1. Comparison of <j)e of PB and CCPB equation with the electroneutral condition. 
4>e of PB equation are in (red) dash-dotted curves and (fe of CCPB equation are in (blue) solid 
curves. The label of curves in each picture depends on the dielectric constant e = 2~^,2~^,2~^. 
I. » 7 e = 0.5e^ and 7 = 0 . II. 7 e = 0.5e and 7 = 0.5. In this computations, three ion species are 
used, one anion with valence — 1 , 0:1 = 1.2 and two cations with valences 1,2 , /3i = /32 = 0.4. 


Table 5.2. The numerical results of </ie(0) of CCPB equation and their limit values t, c, 

c* in El) where (y.\= f3\-\-2f32. /3i[Na"*"] is fixed to 1. e is fixed to 2 


Ve 

P 1 /P 2 

^£( 1 ) 

t 


c 



1 

1.0000 

1.0000 

-0.1124 

-0.1126 

-0.1446 

0 

1/2 

1.0000 

1.0000 

-0.1265 

-0.1265 

-0.1446 


1/3 

1.0000 

1.0000 

-0.1320 

-0.1320 

-0.1446 


1 

0.9679 

1.0000 

-0.1059 

-0.1126 

-0.1446 

0.5e2 

1/2 

0.9581 

1.0000 

-0.1171 

-0.1265 

-0.1446 


1/3 

0.9504 

1.0000 

-0.1206 

-0.1320 

-0.1446 


1 

0.4962 

0.4960 

-0.0299 

-0.0299 

-0.0394 

0.5e 

1/2 

0.4278 

0.4277 

-0.0255 

-0.0255 

-0.0296 


1/3 

0.3853 

0.3853 

-0.0218 

-0.0218 

-0.0242 


From Theorem [O] (ii)-(iv), both t and t — c are decreasing functions to 7 . Surely, 
c can be regarded as a function to 7 . Under some specific conditions, c may become 
a increasing function to 7 (see Remark 11.21 and the graph 1 in each panel in Figure 
ESI). However, it is not clear if the function c has monotonicity generically. Using 
the Newton’s method, we solve the system of equations (ll.lbp and p.l7l) and obtain 
the graph of c and t, respectively. We first consider three ion species with coefficients 
satisfying bi — l,b 2 = 2 and 6i/3i + 62/32 = 01^1 = 1.2. Specific values of /3i and /32 can be 
chosen as follows: 

1. (/3i,/32) = (1.199,0.0005), 

11. (/3i,/32) = (0.002,0.599). 

For each (/ 3 i,/ 32 ), graphs of c and t corresponding to the cases of (ai,ai) = (1,1.2), 
(2,0.6) and (3,0.4) are plotted in Figure lOl respectively. As for Theorem 11.31 lii'l-liiib 
our numerical results indicate that |c( 7 )| < t{j) for all 7 > 0 ; both 0 ( 7 ) and t{j) tend to 
zero as 7 goes to infinity. For each fixed 7 > 0 , the value ^( 7 ) increases but the value 
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c(j) decreases as ai increases. Similar results can also be observed for four ion species 
with coefficients satisfying the following conditions: 

Case 1. aiai=/3i +2/32+ 3/33 = 1.5, (/3i,/32,/33) = (0.25,0.25,0.25), 

01 = 1,2,3,4, i.e., ai = 1.5,0.75,0.5,0.375, 

Case 2. ai + 2a2 = /3i +2/32 = 1.5, (/3i,/32) = (0.75,0.375), 

(ai,(32) = (0.3,0.6),(0.5,0.5),(0.75,0.375), 

Case 3. ai + 2 a 2 = /3i + 2/32 = 1.5, (/ 3 i,^ 2 ) = (0.5,0.5), 

( 01 , 0 : 2 ) = (0.3,0.6),(0.5,0.5),(0.75,0.375), 

Cased. ai+ 2 o 2 = /3i+ 2/32 = 1.5, (/ 3 i,^ 2 ) = (0.3,0.6), 

( 01 , 02 ) = (0.3,0.6),(0.5,0.5),(0.75,0.375). 

The profiles of c and t associated with Case 1-4 are sketched in Figure 15.31 I-IV, 
respectively. As for Figure [521 various a^’s may result in different profiles of function 
c = 0 ( 7 ). However, until now, all our results only show that the function c is of monotone 
increasing or decreasing. This motivates us to see if the function c becomes a non¬ 
monotone function under the other conditions of a^’s and /3j’s. 


I- (+,+) = (1.198,0.001) 



II. (+,+) = (0.002,0.599) 



Fig. 5.2. Comparison of 0 ( 7 ), t(^) with three species; one negative charge, two positive 
charges where ai = 1.2,0.QjOA for 1, 2, 3, respeetively. I. (/3i,/?2) = (1.199,0.0005). II. (/3i,/32) = 
(0.002,0.599). 


As shown in both Figure 15.21 and 15.31 we observe that 0 ( 7 ) converges to zero as 
7 goes to infinity. This is consistent with the results of Theorem 11.31 Moreover, the 
profile of function c can be changed from monotone decreasing to increasing. Such a 
behavior of c and the nonlinearity of equations (I1.16|) and (11.171) let us believe that 
the non-monotone profile of function c may exist. To get the non-monotone profile of 
function c, we consider the following conditions: 

A. 2ai=/3i + 2,a2 + 3/33 = 1.5, (/3i,/32,/33) = (0.9,0.12,0.12), 

B. 2ai=/3i + 2/32 + 3/33 + 4/34 = 1.5, (/3i,/32,/33,/34) = (1.23,0.03,0.03,0.03), 

C. 3ai=/3i + 2/32 + 3^3+4/34 = 1.5, (/3i,/32,/33,^4) = (0.6,0.1,0.1,0.1), 

D. 3ai=/3i + 2/32 + 3/33+4/34 = 1.5, (/3i,/32,/33,/34) = (0.1,0.35,0.1,0.1). 

The non-monotonic profiles of function c with respect to conditions A-D are provided 
in Figure 1231 1-4, respectively. However, the profiles of functions t and t — c are still 
monotonically decreasing. 


6 . Conclusion 

























26 


Boundary layer solutions of Charge Conserving Poisson-Boltzmann equations 


I- (A,ft. A) = (0.25,0.25,0.25) 



7 


II. Qi + 2q2 = ft + 2ft = 1.5, (ft.ft) = (0.75,0.375) 



III. Qi + 2a2 — ft + 2ft — 1.5, (/3i,ft) — (0.5,0.5) 



7 


IV. ui + 202 — fli + 2ft — 1.5, (ft,ft) — (0.3,0.6) 



7 


Fig. 5 . 3 . Comparison 0/0(7), f:(7) with four ion species; one negative charges, three positive 
charges (I), and two negative charge, two positive eharges (II, III, IV). I. ai = 1.5,0.75,0.5,0.375 
for 1, 2, 3, 4, respectively, and (^Si,,02,/Js) = (0.25,0.25,0.25). II. (,0i,/32) = (0.75,0.375). 
III. (^i,/32) = (0.5,0.5). IV. (^i,/32) = (0.3,0.6). (^I.(ai,a2) = (0.3,0.6), 2 .( 01 , 02 ) = (0.5,0.5), 
3 .( 01 , 02 ) = (0.75,0.375) for II, III, IV.) 


For the binary mixture of monovalent anions and cations, although CCPB and 
PB can have very different solutions with different boundary conditions and other con¬ 
straints, the solutions of CCPB equations have very similar asymptotic behavior as 
those of PB equations when the global electroneutrality (11.31) holds (cf. [5D]). 

Situation becomes more complicated in the presence of mixtures of multiple (more 
than three) species with multivalences. In this paper, we again consider the situations 
under global electroneutrality, but the general mixture of multi-species ions. The (more 
rigorous) CCPB shows very different asymptotic behaviors to PB equations under Robin 
type boundary conditions with various coefficients jye’s. 

In particular, the solution i/j of CCPB equation may tend to a constant c at interior 
points, and ±t at boundary points as e goes to zero. As both t and t — c 

are monotone decreasing functions of 7 . Physically, 7 can be regarded as the ratio 
of the Stern-layer width to the Debye screening length. Various conditions can be 
found theoretically and numerically such that the function c of 7 becomes monotone 
decreasing, increasing and non-monotonic. While for PB equation, the solution c/e only 
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Non-monotone profiles of 0 ( 7 ) 



Fig. 5 . 4 . Non-monotone profiles 0/0(7). 


tend to zero at interior points which is independent to 7 . This constitutes one of the 
main differences of PB and CCPB equations. 

This work is one of our first attempts in systematically studying the ionic fluids. 
Much works are needed in the future. In particular, the theoretical justification of the 
interesting behavior of 0 ( 7 ) with respect to 7 under different physical conditions. The 
problems involving multiple spatial dimension domains are for certain to provide more 
interesting phenomena of the solutions and also more technical challenges. Overall, our 
results again demonstrate that the CCPB equation being a more physical and suitable 
model for future applications involving the mixture of multi-species ions. 

Appendix. For the convenience of the readers, we will list out our previous results 
for 2 mono-valence species with charges of opposite signs situations [SD] . 

Considering CCPB equation with Ni= N 2 = 1, ai= 6 i = l, in [30], we had 
established the following results: 

(a) In the electroneutral case (oi =/3i): 

(al) If lim(: 4 ,o^ = 0, the solution (j)^ approaches zero in [—1,1] as e^O. However, 
has slope of order 0{l/%) on the boundary. 

(a2) When ^ > C for some positive constant C independent of e, the solution 
possesses boundary layers with thickness e. 

(b) In the non-electroneutral case (oi ^ fii). 

The solution has boundary layers with thickness and (j)e{x) — (j)e{±l) tends to 
infinity with the leading order term log(e“^) as e^O for x£ (—1,1). The values (^e(±l) 
can be estimated as follows: 

(bl) If ^ < C, </>e(l) and 0g(—1) converge to different finite values as e|0, where C 
is a positive constant independent of e. 

(b2) If lim(:_i,o^ = 00 , both ^e(l) and 0e(—1) diverge to 00 , but |(/)e(l) — 0e(—1)| con¬ 
verges to zero as e^O. 

(c) The difference between the solutions to the CCPB equation m and the PB equa- 
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tion (d can be stated as follows: 

(cl) When ai=/3i, the solution of the CCPB equation (11.21) may converge to the 
solution of the PB equation (11.41) . Namely, in the case of ai =/3i, the solution 
of the CCPB equation dd has the same asymptotic behavior as that of the 
PB equation p.4l) . 

(c2) When ai ^ /3i, the solution of the PB equation (11.41) remain bounded for e > 0. 
However, as ai^Pi, the solution of the CCPB equation dd may tend to 
infinity as e goes to zero (see (b)). This may provide the difference between 
the solutions to the CCPB equation (II.2|) and the PB equation (11.41) . 
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